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We give a concrete experimental scheme for engineering the insulator-superfluid transition of light 
in a one-dimensional (1-D) array of coupled superconducting stripline resonators. In our proposed 
architecture, the on-site interaction and the photon hopping rate can be tuned independently by 
adjusting the transition frequencies of the charge qubits inside the resonators and at the resonator 
junctions, respectively, which permits us to systematically study the quantum phase transition of 
light in a complete parameter space. By combining the techniques of photon-number-dependent 
qubit transition and fast read-out of the qubit state using a separate low-Q resonator mode, the 
statistical property of the excitations in each resonator can be obtained with a high efficiency. An 
analysis of the various decoherence sources and disorders shows that our scheme can serve as a guide 
to coming experiments involving a small number of coupled resonators. 
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I. INTRODUCTION 

In the past two decades, there has been a great inter- 
est in mimicking various quantum many-body phenom- 
ena with artificially engineered structures that permit un- 
precedented experimental control and measurement ac- 
cess [1, 2]. A very successful example in this direction 
is the simulation of Mott insulator-superfluid transition 
with ultracold atoms in optical lattices j3j . Recently, the 
coupled resonator array has been suggested as another 
romising candidate for building a quantum simulator 
5j. Compared to other structures, the coupled res- 
onator array can be used to study the quantum phase 
transition of light and has the striking advantage of full 
addressability of individual sites. 

Since the idea of realizing strongly correlated states 
of light in coupled resonator arrays was first proposed in 
the seminal papers a large amount of work has been 

devoted to a systematic study of the light phase in the 
total parameter space [9l-[l9|. The phase boundary be- 
tween the Mott insulator (MI) phase and the superfluid 
(SF) phase has been obtained using different numerical 
methods including the mean-field app roach 

UM, the 

quantum Monte Carlo simulation [ll|, and the density- 
matrix renormalization-group approach [l^ [isj . Ana- 
lytical and numerical methods were also developed to 
calculate the single-particle excitation spectrum of the 
phase space [lj-[lg|. By considering effects of the dissi- 
pation and driving terms, recent work by several groups 
has promoted the quantum phase transition of light to a 
nonequilibrium case 17-19]. 

Despite the plentiful and substantial achievements in 
the theoretical aspect, there have not yet been any exper- 
imental realizations of coupled resonator arrays. For an 
experimental exploration of the quantum phase transi- 
tion of light using coupled resonator arrays, we should 
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have the abilities of preparing the total system in its 
ground state, tuning the effective photon repulsion and 
the photon hopping rate over a wide range of values, and 
obtaining the accurate statistical property of the excita- 
tions in each resonator. Up to now, previous work has 
just provided some primitive hints toward possible real- 
izations of the model d, [1^. The authors in Refs. 
01 Hi suggested measuring the individual resonator via 
mapping the excitations onto the atomic levels followed 
by state selective resonance fluorescence, but this method 
still suffers from the current lack of high-efficiency photon 
detectors. Moreover, the photon hopping rate is typically 
fixed by the fabrication process and can not be changed 
post-creation, which limits the study of quantum phase 
transition of light in a complete parameter space. 

In this paper, we propose an alternative experimen- 
tal scheme which does not suffer from the above limi- 
tations. For our version of the coupled resonator array, 
each superconducting transmission line resonator (TLR) 
contains a charge qubit as the nonlinear element and 
the adjacent resonators are coupled by another charge 
qubit playing the role of an effective knob for the pho- 
ton hopping rate. The local statistical property of each 
superconducting resonator can be analyzed readily using 
combined techiiiques of photon-number-dependent qubit 
transition [2l|, and fast read-out of the qubit state 
through a separate low-Q resonator mode [23], for which 
the high-efficiency photon detectors are not required. Be- 
cause all the techniques we use have been separately 
demonstrated in the laboratory, our scheme may be im- 
plemented in the near future. 

The paper is organized as follows. In Sec. II we in- 
troduce the physical system considered and derive the 
effective Hamiltonian. In Sec. Ill, a static phase diagram 
in the complete parameter space is obtained using the ex- 
act diagonalization techniques and some related results 
are discussed. Finally, we explore the possibilities of engi- 
neering photonic MI-SF transition dynamically using our 
proposed architecture in Sec. IV. A concrete experimen- 
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FIG. 1. (Color online) (a) Schematic layout of our proposed architecture. A 1-D array of n TLRs R\,R2, ... ,Rn is arranged 
into an annular geometry to satisfy the periodic boundary conditions. Each resonator contains a charge qubit to offer the strong 
nonlinearity. The adjacent resonators Ri and Ri+\ are coupled by another charge qubit d playing the role of an effective knob 
for the photon hopping rate, (b) Electrical circuit representation of the adjacent resonators and their junction. The high-Q 
half- wave mode of Ri with resonance frequency w is coupled to transmon i, and RiS low-Q full- wave mode with resonance 
frequency 2'w is strongly coupled to a measurement line fabricated at the resonator center. The state of transmon i can be 
measured by applying a microwave field (7™ at the input port of the measurement line. Another transmon d is injected into 
the resonator junction , which is dispersively coupled to Ri and -Ri-i-i. Microwave pulse Uf applied to the gate serves for driving 
the qubit transition of transmon i. The transition frequency of each transmon can be tuned via the corresponding applied 
magnetic flux $. 



tal procedure is presented including how to initialize the 
system, tune the system's parameters, and measure the 
individual resonators. The various decoherence sources 
and disorders are also analyzed. 



II. THE PHYSICAL SYSTEM AND EFFECTIVE 
HAMILTONIAN 

The system we consider is schematically depicted in 
Fig. 1(a). A 1-D array of n TLRs i?i,i?2, ■•■ ^Rn is ar- 
ranged into an annular geometry to satisfy the periodic 
boundary conditions. Each resonator contains a charge 
qubit (represented by the red dot) to offer the strong 
nonlinearity. The adjacent resonators Ri and Ri+i are 
coupled by another charge qubit Ci (represented by the 
blue cross) playing the role of an effective knob for the 
photon hopping rate. The charge qubit used in our model 
is transmon, a modified version of the Cooper pair box 
proposed by Koch et al.^^. Its unique feature is the 
large shunt capacitor between the superconducting is- 
lands, which makes the transmon have a longer deco- 
herence time than the ordinary charge qubits. High-Q 
resonators are advantageous for the simulation of MI- 
SF transition but are adverse to the measurement of in- 
dividual sites. To solve this problem, we can use the 
technique of engineering two modes of a resonator with 
different quality factors, which has been demonstrated 
experimentally in [2^. As shown in Fig. 1(b), the high-Q 
half-wave mode of Ri with resonance frequency w is cou- 
pled to transmon i, and RiS low-Q full- wave mode with 
resonance frequency 2w is strongly coupled to a measure- 
ment line fabricated at the resonator center. The state 
of transmon i can be measured by applying a microwave 
field [/™ at the input port of the measurement line. An- 
other transmon Ci is injected into the resonator junction. 



which is dispersively coupled to Ri and Ri+i- Note that 
the usual shunt capacitor between the superconducting 
islands of Ci is replaced with capacitors to the ground 
planes to suppress direct coupling between Ri and Ri+i 
[23 |. Microwave pulse Uf applied to the gate serves for 
driving the qubit transition of transmon i. The transi- 
tion frequency of each transmon can be tuned via the 
corresponding applied magnetic flux $. 

In the following, we will derive the effective Hamil- 
tonian of our proposed architecture. Let us denote the 
lowest two eigenstates of transmon i with \g)i and |e)i, 
which are separated by energy e and coupled to the half- 
wave mode of Ri with qubit-resonator coupling strength 
g. CiS lowest two eigenstates \g'^)i and \e'^)i, separated 
by energy Ec, are simultaneously coupled to Ri and Ri-^-i 
with coupling strength gc- In this paper, our analysis is 
restricted to the case oi ec — w ^ gc, i- e. , Ci is disper- 
sively coupled to its neighbor resonators. The Hamilto- 
nian for the total system can be written as (assuming 



H = Hi+H2, 



(e| + wajfli + 



gia+ai 



4)], 



H2 = ^[ee|e=)..(e=| + gc{<J+a, + a+a,+^ + H.c.)], (1) 
1=1 

where Hi is the sum of local Jaynes-Cummings Hamil- 
tonians with resonator index i, photon creation (annihi- 
lation) operator a\ (ai) and qubit raising (lowering) op- 
erator af' (cr~) for transmon i; H2 describes the sum of 
interactions between Ci{i — 1, 2, . . . , n) and their neigh- 
bor resonators with ct^ i'^ci) being C^'s qubit raising 
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(lowering) operator. Considering the annular geometry 
of the coupled resonator array, we have a„+i — Oi and 
al^_^_i = a\, which offer the periodic boundary conditions 
for our architecture. 

In an interaction picture with respect to Hq = 

Sr=i[^(l'3)"('5| + "'l^i) + '^c\e'')u{e''\], the system Hamil- 
tonian reads 

n 

Fr = Y.^A\e)u{e\+9{a+a, + a-al)], 

i=l 



where A' = A + and k(Ac) = is the Ac- 
dependent photon hopping rate. 

The Hamiltonian H^'^p^' in Eq. (5) actually provides 
an effective (A, Ac)-dependent on-site repulsion for exci- 
tations. H^'^P'^^ can be diagonalized in a basis of mixed 
photonic and atomic excitations, called polaritons[^ 0]. 
Let \n,g) e)) represent a resonator that contains n 
photons and a single qubit in the ground (excited) state. 
Then the polariton states of Ri , labeled by the polariton 
number n and upper or lower branch index cr = ±, can 
be given by 

|n, +)i = sin6'„|n - 1, e)i + cos 9n\n, g),, 



n 

HT' = ^[5c^+(a.+a,+i)e'^=*-H5c^;,(aI + aI+i)e-'^^*], 
i=i 

(2) 

where A = e — w, and Ac = ec — w. With the choice 
of Ac 3> gci the real energy exchanges between Ci{i = 
1, 2, . . . , n) and their neighbor resonators are largely sup- 
pressed. In this case, we can use the time-averaging 
method in Ref. [1^ and neglect the effect of rapidly os- 
cillating terms. Then can be approximated by 

" 1 

■ ■ 1 

n 2 

+al^{a\^^a, -f H.c.) + [a+cJ-^^^ + H.c.)], (3) 

where al^ = \e'^)ii{e'^\ — \g'^)ii{g'^\- The four terms 
of H2 describe the ac Stark shifts, Lamb shifts, 
transmon-intermediated photon hoppings, and photon- 
intermediated dipole couplings between the transnions 
situated at the neighbored junctions, respectively. If 
Ci{i = l,2,...,n) are prepared in \g'')i\g'')2 ■ ■ ■\9'')n, 
then they will always stay in their ground states. In 
this situation, iJg can be simplified by 

= „(.9l . . .2 (ff1l(/|i/;i5^)l|.g^)2 . . . 

= Yli-'^ial+iai + aU^+l) - -^4ai]- (4) 

z— 1 ^ ^ 

By moving to a second interaction picture with respect to 
Sr=i + kliiieDj the total system Hamiltonian 

j"* + i?2 yields 

n 

jjhop ^ _^(Ac)(al+ia, + ata,+i), 

n 

jjrepui ^ Y,[A\e}u{e\+g{ata. + a-al)], (5) 

i=l 



|n, -)i = cos6'„|n - 1, e)j - sm6n\n, g)^, (6) 



with the mixing angle tan6'„ ~ — — ^y=^ • The 

corresponding eigenvalues are 

The zero-polariton state |0, — )i = \0,g)i is a special 
case with Eq = 0. Obviously, these polariton states 
are also eigenstates of the polariton number operator 
Ni — a|ai -I- |e)ii(e| with eigenvalue n. For the case of 
k(Ac) not much larger than g, if we create only the lower 
polariton states in the resonator array initially, then the 
upper polariton states will never be created because the 
interconversion between the two polariton branches can 
be neglected. The effective on-site repulsion Uef / for po- 
laritons results from the anharmonicity in the spectrum 
of iJ'"'^P"'^ which is dependent on the number of polari- 
tons in the resonator. For polariton number n ^ I, the 
effective repulsion Ueff{l) can be identified by calculat- 
ing the energy cost to inject a second polariton into the 
resonator, 

C/c// (1) = - 2E^ 

= -^ + V^A'VV-y(^)2 + 2g2. (8) 

Ueff{l) can be tuned easily by changing the detuning 

a' = A -I- If a' > and a' > 5, C/e//(l) vanishes; 

if a' < and —A' g, then Ueff{l) ~ —A' is a large 
quantity. 

III. STATIC PHASE DIAGRAM IN THE (A, A^) 
PARAMETER SPACE 

In Sec. II, we have derived the effective Hamiltonian 
^e// |.]-^g proposed architecture, which has two in 
situ tunable parameters A and Ac. In different regimes 
of the (A, Ac) plane, the system can show distinct char- 
acteristics. For simplicity, our analysis is restricted to 
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the case of the resonator array contains on average one 
polariton per resonator. 

o 2 2 

If A + = and ^ ^ g, the system Hamiltonian 
is dominated by the resonant Jaynes-Cummings interac- 
tion, i.e. H'^f^ « X]"=i + O'^^ local res- 
onator has a polariton in it, the strong qubit-resonator 
interaction will shift the frequency of the resonator mode 
and prevent a second polariton from entering it [26l |. 
This anharmonicity in the spectrum leads to an effec- 
tive polariton-polariton repulsion Ueff{l) = {2~y/2)g ^ 

2 

k(Ac) = In this case, the ground state of the system 
is approximately 



MI 



(9) 



which can be seen as the MI state of polaritons. 
If A = A + ^ ^ g, we have \n, —)i « 



~\n,g)i 

and C/e//(l) ~ 0. When k(Ac) = ^ is nonvanish- 
ing, the photon hopping terms will dominate the system 

Hamiltonian, i.e-H^^f « H'^^p = - Er=i ^(al+i^* + 
a|ai+i). In this situation, H'^f-^ can be diagonal- 
ized through the Fourier transform. Introducing bk = 
E"=i;^e'^aj (fc = 0, 1,2, . . . ,n - 1), we have 

[bk,bl] = 1 and H'^ff can be rewritten as H^^f « H'"'p = 

"^S^ X]fe=o '^os ^^6|,6fe. Obviously, the system ground 
state is given by 



^)SF = ^{blrivac) 



1 . 1 



— i—yalr\vac) 



/ J I J 



(10) 



and the corresponding ground state energy is — 2n-^, 
where \vac) is the vacuum state of the resonator array. 
With n photons delocalizing throughout the whole res- 
onator array, |^')s_f is the SF phase state of light. 

By calculating the system ground states correspond- 
ing to different values of A and Ac, we can obtain a 
static phase diagram in the (A, Ac) plane and observe 
the phase boundary between the MI phase and the SF 
phase. In this section, we will show that, even with a very 
small resonator array, the phase diagram can exhibit the 
main features of MI-SF transition, which is very advanta- 
geous for the potential experimental realization. To plot 
the phase diagram, we must choose an order parameter 
to differentiate between insulatorlike and superfluidlike 
states. For small size resonator arrays, a good order pa- 
rameter is the variance of the polariton number in a given 
site [7, 8], which is defined as 



Eq. (9), var{N.C) = 0; in the case of n = 3, the SF 
state |5')5F in Eq. (10) results in a nonvanishing variance 
var{N,) « 0.6668. 

The MI-SF phase transition results from the interplay 
between the photon hopping and the on-site repulsive in- 
teraction. As shown in Fig. 2(a), if we choose the system 
parameters g = ^c, A e [— 5g, IO5] and Ac € [IO5, 100(/], 
the ratio of the photon hopping rate to the on-site repul- 
sion K/Ueff{l) can be tuned from smaller than 10^^ to 
larger than 10, which permits us to systematically study 
the phase transition in the (A, Ac) plane. Note that our 
choice of Ac satisfies the dispersive coupling condition 
Ac ^ (?c, which has been used in deriving the effective 
Hamiltonian Eq. (5). 

For a 3-site resonator array containing on average 
one polariton per resonator, the system basis consists 
of state vectors of the form \(j))k = (^i^iln^ , s'^)i, where 

e {0,1,2,3}, s'^ e {g,e} and fe(</>|(E ti = 3. 

It is easy to check that there are a total of 38 state vectors 
in this form. To plot the phase diagram in (A, Ac) plane, 
we must calculate the order parameter var{Ni) of the sys- 
tem ground state as a function of A and Ac as follows. 
First: obtain the matrix elements of iJ^-^^ in the system 
basis vectors {k — 1,2,..., 38). Second: diagonalize 
this 38 X 38 matrix and identify the eigenvector \<j>)g cor- 
responding to the lowest eigenvalue. Third: compute the 
order parameter var{Ni) =g {(j)\Nf\(j))g - {g{(l}\N,\(j))g)^. 
The obtained values of var{Ni) as a function of A and 
Ac are plotted in Fig. 2(b). It is shown that, although 
the size of our resonator array is very small, the phase 
transition takes place over a narrow variation of the 
parameters. Similar to the Bose-Hubbard model, the 
value of var(Ni) is mainly determined by the ratio of 
the photon hopping rate to the on-site repulsion [131 ■ 
For our 3-site resonator array, the phase boundary be- 
tween MI and SF can be approximated by the contour 
line K/Ueff{l) = 0.28. By fixing A at a certain value and 
varying Ac, or fixing Ac and varying A, we are all able 
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var{Ni) ^< Nf > ^ < N, >^ 



(11) 



In the MI phase, the number of polaritons per resonator 
is well defined and has zero variance. However, in the 
SF phase, the polariton number fluctuates and thus the 
variance has a nonzero value. Obviously, for |5')Af/ in 



FIG. 2. (Color online) (a) Contour lines of logj^Q(K/C/e//(l)) 
in the (A, Ac) parameter space. The system parameters we 
choose are g — g^, A £ [—5g,10g], and Ac £ [lOg, lOOy]. 
(b) Order parameter var{Ni) of the system ground state as a 
function of A and Ac for a 3-site resonator array. The phase 
boundary between MI and SF can be approximated by the 
contour line n/UejfiX) ~ 0.28(black solid line). 
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to tune the system from MI phase to SF phase. How- 
ever, the latter method is more effective because var{Ni) 
is more sensitive to A. 

The phase diagram presented above is obtained us- 
ing the effective Hamiltonian H'^^^ , which neglects the 
real energy exchanges between Ci and their neighbor res- 
onators. Because the off-resonant transitions of Ci occur 
with the probability ^ (-2^)^, the Hamiltonian H'^^^ be- 
comes accurate when increases. In Fig. 3, we compare 
the order parameters of the approximate and exact sys- 
tem ground states for a 3-site resonator array, obtained 
from numerical diagonalization of the effective Hamilto- 
nian H'^f-f and the full Hamiltonian H in Eq. (1), respec- 
tively. The parameters we choose are g — gc, Ac — lOgc, 
and A g [—5(7, lOg]. The good agreement between these 
results indicates that Ac = lOgc is large enough to guar- 
antee the accuracy of H'^-f^ . 



as the number of sites is increased. 




FIG. 3. (Color online) Order parameters of the approximate 
and exact system ground states for a 3-site resonator ar- 
ray, obtained from numerical diagonalization of the effective 
Hamiltonian H"^^ (black solid line) and the full Hamiltonian 
H in Eq. (1) (red crosses), respectively. The parameters we 
choose are g = gc, Ac — Wgc, and A G [— 5g, Wg]. 

The increase in the number of sites will lead to a 
sharper phase transition (l^ . Suffering from the expo- 
nentially growing Hilbert space with the resonator array 
size, the exact matrix diagonalization method we use here 
can only obtain the solutions for very small resonator ar- 
rays. To give some primitive hints for the relation be- 
tween the sharpness of the phase transition and the res- 
onator array size, we compare the order parameters of 
the system ground states for 2-site, 3-site and 4-site res- 
onator arrays in Fig. 4. The parameters we choose are 
g = gc, Ac = lOgc, and A e [~5g,10g]. For simplic- 
ity, the Hamiltonian we use is the effective Hamiltonian 
H'^f^. It is shown that, the MI-SF transition takes place 
at the almost same value of A for different resonator ar- 
rays, and we may observe a more abrupt phase transition 
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FIG. 4. (Color online) Order parameters of the system ground 
states for 2-site (blue dashdotted line), 3-site (red dashed line) 
and 4-site (black solid line) resonator arrays. The parameters 
we choose are g — gc, Ac — Wgc, and A £ [~5g, Wg]- 

As presented avove, with a small size system of our 
proposed architecture, we may observe the main charac- 
ters of MI-SF transition experimentally. In next section, 
we will discuss the related experimental issues. 



IV. DYNAMICAL OBSERVATION OF THE 
PHOTONIC MI-SF TRANSITION 

In the following, we give a concrete experimental pro- 
cedure to engineer photonic MI-SF transition dynami- 
cally using our proposed architecture. We will describe 
in detail the total manipulation process (initialization, 
evolution and measurement) and analyze the various ex- 
perimental imperfections. 

A simple scheme to initialize the system is the res- 
onant pumping approach [Sj, 20^. We start in the MI 
regime with the total system in its absolute ground state 
^?=i\0,g}i- Then, by applying a global external mi- 
crowave TT-pulse, we are able to drive the filling factor 
of each resonator from zero to one and prepare the sys- 



tem in the required state 



1, 



With the current 



circuit QED experimental techniques, this preparation 
method has the high fidelity larger than 99%. 

Now, we illustrate how to measure the variance of the 
polariton number in the resonator Ri. Our procedure 
utilizes two harmonic modes of Ri which are engineered 
to have very different quality factors As shown in 

Fig. 1(b), Ri's full- wave mode with resonance frequency 
2w has an electric field antinode at the resonator cen- 
ter, and is hence strongly coupled to the measurement 
line. By choosing big coupling capacitances for the mea- 
surement line, the full-wave mode has a low quality fac- 
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tor. Conversely, Ri's half- wave mode with resonance fre- 
quency w has an electric field node at the resonator center 
and couple weakly to the measurement line, leaving the 
quality factor of this mode limited only by internal losses. 
To obtain var{Ni), we only need to get the probability 
distribution (I — 0,1,2, ... , n) of the polariton num- 
ber, withp, the probability of finding exactly I polaritons 
in Ri. In the experiment, one can obtain the value of 
as follows: (1) First, switch off the effective polaritonic 
hopping instantaneously to isolate the system state from 
further evolution (this can be achieved by tuning the de- 
tuning Ac to a very large value instantaneously). (2) 
Adjust the detuning A adiabatically such that A/g ~ 5. 
At this detuning, the polaritons are transferred into mi- 
crowave photons, i.e.\n,—)i « —\n,g)i, and the qubit 
frequency is strongly shifted depending on the number of 
photons in the resonator (2l| . (3) Then, we drive trans- 

2 

mon i at the frequency (e -I- 2/^), thus to selectively 
populate the qubit into the excited state \e)i if there are 
I photons in the resonator. (4) Next, tune the qubit tran- 
sition frequency e such that transmon i is decoupled from 
i?i's half- wave mode but dispersively coupled to Ri's full- 
wave mode. Then, the state of transmon i can be mea- 
sured fast by applying a microwave field U^^ of frequency 
2w at the input port of the measurement line [23| . (5) 
Repeat steps (l)-(4) for a large number of times M, and 
obtain the number of times Mi in which transmon i is 
excited into \e)i. Provided that we have always prepared 
the same system state before every implementation of 
steps (l)-(4), p, can be obtained as p, « Finally, we 
can get the order parameter as 

n n 

var(K{)=Y.^\~{Y.^p,f. (12) 
/=o /=o 

For our annular architecture, performing the transmis- 
sion measurement of individual resonators requires fabri- 
cating bridges on the circuit. To avoid this experimental 
difficulty, we can improve the measurement procedure by 
measuring reflection amplitudes instead of transmission 
amplitudes [22i]. 

Using the proposed architecture, we are able to observe 
the dynamical quantum phase transition of microwave 
photons by adiabatically changing the photon hopping 
to the on-site repulsion ratio. Similar approach has been 
adopted in the optical lattice experiment Now, we 
analyze the feasibility of this procedure by some rough 
calculations based on the practical experimental parame- 
ters. It has been shown that, to make the system always 
remain in the many-body ground state of the varying 
Hamiltonian, the timescale of tuning the system from the 
MI state to the SF state or vice versa should be compara- 
ble to the photon hopping time [1]. If we perform the 
experiment by fixing Ac at a certain value and varying 
A, the photon hopping time ^ = Coupling strength 
cjc = 27r X 200 MHz has been realized experimentally 
in Ref. For Ac € [lOgc, lOOgJ, the photon hopping 
time can be calculated as i = 8 ~ 80 ns. On the other 



hand, an estimate of the polariton lifetime Tp = 2 fis 
is reasonable for the practical situation |23j]. It is clear 
that the timescale of the many-body dynamics is much 
shorter than the decay time of the mixed photon-qubit 
excitations, so our proposal is feasible with present-day 
technology. 

In the practical experimental realizations, some system 
disorders are unavoidable because of lithography errors. 
The difference of the qubit transition frequencies can be 
compensated by controlling the local magnetic fiuxes $. 
The variation of the resonator frequencies is relatively dif- 
ficult to control. To solve this problem, we can compen- 
sate this variation by inserting a superconducting quan- 
tum interference device (SQUID) into the central con- 
ductor of each resonator to make its frequency tunable 
(28j . However, such an injected SQUID will introduce 
anharmonicity into the resonator energy levels. As a re- 
sult, the effective on-site repulsion is composed of two 
parts: one part comes from the local Jaynes-Cummings 
interaction, the other part comes from the introduced 
resonator anharmonicity. Denoting the on-site repulsion 
induced by the resonator anharmonicity with C^e//' 
critical point for the phase transition will be changed 
from 77-^ ^ 0.28 to - — -^^—r- ~ 0.28. To observe 

the MI-SF transition, U^j^j: has to fulfill the condition 

^'eff 0%' With the current fabrication techniques, 
the coupling strengths typically have larger relative fluc- 
tuations than the resonator frequencies. The variation 
in g is harmless because the on-site repulsion [/e//(l) 
is a complicated function of g and A, and the errors in 
C/e//(l) induced by the deviation of g can be compen- 
sated by individually adjusting the detuning A. How- 
ever, if the coupling strength gc fluctuates from sample 
to sample, the ac Stark shifts of the resonator frequencies 
in Eq. (3) can not be removed by going to the interaction 
picture. These terms are equivalent to additional param- 
eter spreads in the resonator frequencies. Consequently, 
both the deviations in the resonator frequencies and the 
coupling strengths can be compensated by making the 
resonators tunable except it introduces additional part 
of effective on-site repulsion f^e//- 

In conclusion, we propose to observe the dynamical 
quantum phase transition of light using a 1-D array of 
coupled superconducting TLRs. The unique feature of 
our architecture is the good tunability of the photon 
hopping rate and the effective on-site repulsion. The lo- 
cal statistical property of each superconducting resonator 
can be analyzed readily using existing microwave tech- 
niques. With the recent progress in the multi-resonator 
experiments [2^, our proposal may serve as a guide to 
coming experiments of quantum phase transition based 
on a small-scale resonator array. 
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